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We investigate two simple prescriptions to account for the Pauli principle in a three-body cluster 
model employing a new method based on an adiabatic hyperspherical expansion to solve the Faddeev 
equations in coordinate space. The resulting wave functions are computed and compared. They 
\& ■ are furthermore tested on halo nuclei by calculations of momentum distributions and invariant 

mass spectra arising after fragmentation of fast 6 He and 11 Li in collisions with light targets. The 
prescriptions are very accurate and the available measured quantities are remarkably well reproduced 
when final state interactions are included. 
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I. INTRODUCTION 



Three-body models are useful to describe the relative wave function of three particles when their intrinsic structure 
£C) , remains unchanged Q-j^]. The intrinsic and the relative degrees of freedom are then assumed to be completely 
decoupled. However, when the particles themselves are composite structures (clusters) and contain identical fermions 
this assumption violates the Pauli principle related to fermions distributed in different clusters. This violation typically 
leads to unphysical deeply bound two-cluster states. In the two-body problem this violation can be easily cured by 
removing this unphysical state from the active space available for the system. 

In the three-cluster problem most of the usual approaches to include the Pauli principle basically fall into two 
categories: i) to project out the undesired overlap of the three-body wave function with the Pauli forbidden two- 
body state [0-p[ and ii) to modify the two-body interaction in such a way that the unphysical bound cluster-cluster 
state is avoided, provided the low energy scattering properties or the phase shifts remain unchanged. The 
prescriptions basically amount to either excluding occupation of Pauli forbidden two-body states in the three-body 
wave function or excluding the apriori appearance of such undesired states. The two prescriptions give somewhat 
different results and it is not yet clear which one is preferable. The projection technique, although seems to be more 
straightforward, is in practice technically more difficult. 

The problem of treating the intercluster Pauli principle again attracted attention in investigations of the structure 
of the three-body halo nuclei where most of the information is obtained by studying momentum distributions of 
fragments of these nuclei in fast collisions with light targets. Originally these momentum distributions were expected 
to reveal the structure of the initial three-body wave function ]10| . This is still true although more elaborate analyses 
are required. The correct account for the two-body interaction between fragments in the final state turned out to be 
essential when low lying resonances are present Jll| ) p2[ . 

The purpose of the present paper is twofold. First we describe how the projecting out technique can be easily 
incorporated in the adiabatic hyperspherical expansion of the Faddeev equations in coordinate space |13j1 . This 
method to solve the Faddeev equations has already been used in several investigations of halo nuclei f| |l4||'5[ . 

Secondly we investigate how the differences in the prescriptions are reflected in the momentum distributions of frag- 
ments from the break-up reactions of three-body halos. This question is important for an unambiguous interpretation 
of the wealth of experimental data on three-body halos like 6 He ( 4 He+n+n) and 11 Li ( 9 Li+n+n). 

After the introduction we shall in section 2 describe the method and compare potentials and wave functions for 
two different prescriptions. In section 3 we use these wave functions in computations of momentum distributions in 
nuclear break-up reactions of three-body halo nuclei. Finally in section 4 we give a brief summary and the conclusions. 

II. THREE-CLUSTER PROBLEM WITH PAULI FORBIDDEN TWO-CLUSTER STATES 

Adiabatic hyperspherical expansion of the Faddeev equations. We shall use the hyperspherical coordinates (p, a, 
Q x , fly) defined in refs. |^||- The volume element is given by p 5 dQdp, where dfl — sin 2 a cos 2 adad£l x dVL y . The 
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total wave function ^ jm of the three-body system (with total spin J and projection M) is written as a sum of three 
components, which in turn for each hyperradius p are expanded in a complete set of generalized angular functions 
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where p 5 / 2 is related to the volume element. 

These wave functions satisfy the angular part of the three Faddecv equations: 



2mp 2 " + + <± W + <i, n ^- 2 mp 2 " W " ' 



where {«, J, A:} is a cyclic permutation of {1, 2, 3}, m is an arbitrary normalization mass, and A 2 is the p-independent 
part of the kinetic energy operator. The analytic expressions for A 2 and the kinetic energy operator can for instance 
be found in ||. 

The radial expansion coefficients f n (p) are obtained from a coupled set of "radial" differential equations ||, i.e. 

d 2 2mE 1 / . . n 15\\ . . . 

— " — — + -o A„(p) - Qnn + T Mp) = 



dp 



where the functions P and Q arc defined as angular integrals 
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E (2Pnn'-^ + Qn„')/n'(p) , (3) 



p nn ,(p) = E / dfi*<?>,fi) jU#(p,fi) 5 (4) 



Q nn /( P ) = E / ^^'(p.^^^Cp.n) ■ (5) 

»,i=i 



The radial equations in eq.(|3j) reveal that h (A n (p) — Q„„ + j)/(2mp ) is the diagonal part of the effective radial 
potential. Its behavior is decisive for the properties of the three-body system. The A-spectrum at both p = and 
p = oo is identical to the hyperspherical spectrum. In addition, for every bound two-body state, there exists one 
A- value which bends over and diverges parabolically as function of p for p ~ > oo, see Such a level corresponds 
at large distances to the three-body structure where the two-body subsystem is in the corresponding bound state, 
whereas the third particle is far away. 

Treatment of the Pauli principle. The first prescription is a two-body potential without Pauli forbidden states 
but with the same low energy properties as the original deep two-body potential with the forbidden state. This 
corresponds to the use of an additional Pauli repulsion. We have discarded the apparently rigorous procedure of 
constructing the phase equivalent potential from the original two-body potential. The reasons are that the original 
interaction in any case appears in the form of an effective (mean) potential adjusted to reproduce specific measured 
quantities. Furthermore, the binding energy of the three-body system must be accurately reproduced to provide the 
correct structure. An additional fine tuning is therefore almost inevitable. The strictly phase equivalent potential is 
substantially more difficult to obtain and use. Under these circumstances we prefer to take the pragmatic approach 
and construct phenomenological potentials adjusted to have specifically chosen crucial properties. 

The second prescription is a unique feature of the adiabatical expansion method. Since in this method each of 
the two-body bound states gives rise to a separate diverging A-value, the prescription is simply to omit the A-values 
corresponding to the Pauli forbidden two-body bound states from the expansion (|l|). Then all the Pauli forbidden 
states are automatically excluded in the three-body wave functions. 

III. TWO-BODY POTENTIALS 

We shall consider the two halo nuclei 6 He and 11 Li, both approximated as three-body systems consisting of a core 
( 4 He and 9 Li) surrounded by two valence neutrons. This approximation works remarkably well for these systems and 
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furthermore due to relatively weak binding of the systems only the low energy parameters of the two-body potentials 
are essential for the structure of the systems 0]. As the core is inert in our approximation it is sufficient for our 
purpose to assume zero spin of the core for 11 Li. 

In both cases the lowest neutron s-states in the core are filled. Therefore to conform with the Pauli principle the 
halo neutron should be prevented from occupying the lowest neutron-core s-state. We achieve this goal by use of 
two different neutron-core potentials with the same low energy properties. One of the potentials will have a Pauli 
forbidden bound state the other will not. This Pauli forbidden state will be consequently removed from the active 
space. 

Neutron-neutron potential. As indicated in Q and proved by our test runs the particular radial shape of the 
n — n interaction is not important for 6 He and n Li ground state properties as l ong as the low energy n — n scattering 
parameters are correct. We therefore use a simple potential similar to the one in |15| which reproduces the experimental 
s- and p-wave scattering lengths and effective ranges. It contains central, spin-orbit (L • S), tensor (<Si 2 ) and spin-spin 
(s! • s 2 ) interactions and is explicitly given as 

V nn {r) = 37.05 exp(-(r/1.31) 2 ) - 7.38 cxp(-(r/1.84) 2 ) 
-23.77 exp(-(r/1.45) 2 )L • S + 7.16 exp(-(r/2.43) 2 )S , i 2 
+ (49.40 exp(-(r/1.31) 2 ) + 29.53 exp(-(r/1.84) 2 )) Si • s 2 , (6) 

where the strengths are in MeV and ranges in fm. Its scattering lengths!] a and effective ranges r e are (in fm) a( 1 So)= 
18.45, r e ( 1 S )= 2.83, a( 3 P ) = 3.38, r e ( 3 P )= 1.10, a( 3 Pi)= -2.02, r e (^Pi)= -2.94, a( 3 P 2 )= 0.31, r e ( 3 P 2 )= 18.73. 

Neutron- 41 He potential. For 6 He the low energy properties of the neutron-core subsystem in the s- and p- waves are 
rather well known and reflected in an s-wave scattering length of a s — — 3.07± 0.02 fm and two resonance energies E 
and widths V of E(p 3/2 ) = 0.77 MeV, T(p 3/2 ) = 0.64 MeV and E(p 1/2 ) w 1.97 MeV, T(p 1/2 ) w 5.22 MeV, see @. 
The P1/2-P3/2 splitting demands a spin-orbit force or equivalently (with the same range of the interactions) different 
strengths in the P1/2 and P3/2 channels. The negative s-wave scattering length can be reproduced by a repulsive 
potential or by an attractive potential with a bound state which in this case is the Pauli forbidden state occupied by 
the core neutrons. 

We introduce central 1-dependent and spin-orbit components in the potential with a gaussian shape S'exp(— r 2 /b 2 ). 
The range was chosen to be b = 2.33 fm for all components (similar to except for the repulsive s-wave where 
the range b was changed to 3.34 fm in order to reproduce the same effective range as for attractive potential. The 
strength parameters are then defined by the fit to the specified scattering length and positions of the resonances. 

The p-wave strengths are S(pi/ 2 ) = —48.675 MeV and S(p 3 / 2 ) = —53.175 MeV (or, equivalently, the central 
p-wave strength of -51.675 Mev and the spin-orbit strength of -3.00 MeV) which provide the p-wave resonances 
E(p 1/2 ) = 1.94 MeV, T{p 1/2 ) = 4.0 MeV and E(p 3/2 ) = 0.77 MeV, T(p 3/2 ) = 0.73 MeV. These parameters are the 
same for both repulsive and attractive s-wave potential. 

The s-wave strength for the attractive potential (further referred to as "attractive") is <S'(s 1 / 2 ) = —39.2 MeV 
(r e = 1.41 fm and a = —3.07 fm, one Pauli forbidden bound state). Without additional fit these s- and p-wave potentials 
provide the binding energy B( 6 He)=1.0 MeV and r.m.s. radius R( 6 He)=2.45 fm, which is close to experimental data 
B( 6 He)=0.97± 0.04 MeV and R( 6 He)= 2.57 ± 0.10 fm. 

The repulsive s-wave potential which reproduces the experimental scattering length and the same effective range 
as the attractive potential slightly underbindes 6 He (by approximately 200 keV). Such underbinding for repulsive 
potentials is not unusual for 6 He. To alleviate this problem people normally increase the range of all potentials by a 
few per cent However this procedure shifts the positions of the p-resonances from their experimental values. As 
we shall see below these positions are extremely important for the momentum distributions. We therefore leave the 
p- waves unchanged and instead reduce slightly the repulsion in the s-wave. 

The resulting potential is further referred to as "repulsive" and has the s-wave range and strength b = 3.34 fm, 
S{s 1/2 ) = 9.70 MeV, (a = -2.58 fm, r e =0.67 fm) which leads to the binding energy B( 6 He)=1.0 MeV and r.m.s. 
radius R( 6 He)=2.50 fm. 

Neutron- 9 Li potential. For u Li the low energy neutron-core data are less known, although evidence is accumulating 
for a low lying virtual s-state at E(s) w 0.15 ± 0.15 MeV and the lowest p-resonance at E(p) w 0.6 ± 0.2 MeV, see 
Jl7| p0[ . The Pauli forbidden states are in this case both the lowest s-state and the p 3 / 2 -state. For simplicity the latter 
is in the calculations placed at high positive energy and thereby removed from the active space by a large inverse 
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spin-orbit potential Jig ] . For the s- waves we use a shallow potential without bound state and a deep potential with a 
Pauli forbidden state. 

The radial shapes of the neutron-core interactions are also assumed to be gaussians, i.e. 5exp(— r 2 /b 2 ) with 
b = 2.55 fm except for the deep s-wave potential with a bound state where we use b = 1.49 fm to maintain the 
same effective range as for the shallow potential. 

We adopt further the usual assumption that the neutron-core interactions do not depend on the spin of the 9 Li-core. 
A more realistic study of 11 Li properties taking into account the spin of 9 Li has been made in |0,ff5l. 

Again the spin-orbit neutron-core potential effectively only gives different strengths for the two different p-waves. 
With the choice of range for the radial potentials there are only three strength parameters left each related to a 
resonance, a virtual state or a scattering length. The spin-orbit force is used to remove the Pauli forbidden p 3 / 2 -state 
from the active space. The two remaining strength parameters then determine the s- and pi/2-state as well as the 
binding energy of n Li. One of them must be used to fine tune the binding energy B of the total three-body system 
B( n Li)= 0.295 ± 0.035 MeV. With correct binding energy the root mean square radius is then always reproduced 
within the experimental uncertainty, which is R( 11 Li)= 3.1 ± 0.3 fm. 

The potential with a bound Pauli forbidden s-state (further referred to as "deep") has the parameters b = 1.49 fm 
for the s-wave and b = 2.55 fm for p-waves, S(s 1 / 2 ) = —176.608 MeV (one bound state, a=8.738 fm, vir- 
tual level at 0.20 MeV), S(p 3/2 ) = 9.55 MeV (E(p 3/2 ) is high and uninteresting) and S(p 1/2 ) = -38.34 MeV 
{E(p 1 / 2 ) = 0.77 MeV, F(p 1 / 2 ) = 0.89 MeV). The binding energy and root mean square radius is then computed 
to be B( n Li)=0.30 MeV and R( n Li)=3.34 fm. 

The shallow potential without bound states but with the same low energy properties as the deep potential overbinds 
n Li by some 180 keV. In contrast to 6 He the dominating component is now the s-wave in the neutron-core subsystem. 
Therefore to fine tune the binding energy we modify now the p-wave potential namely the position of the pi/ 2 resonance 
while keeping the other properties unmodified. 

The fine tuned potentials without a bound state (further referred to as "shallow") has the parameters b = 2.55 fm, 
S(si/ 2 ) = -7.14 MeV (no bound state, a=8.738 fm, virtual level at 0.20 MeV), S(p 3/2 ) = 9.55 MeV (E(p 3/2 ) is high) 
and S{pi/ 2 ) = -35.45 MeV {E{p 1/2 ) = 1.7 MeV, T(p 1/2 ) = 4.2 MeV). The binding energy and root mean square 
radius is now computed to be B( n Li)= 0.30 MeV and R( n Li)= 3.35 fm. 
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FIG. 1. Angular eigenvalue spectra \ n as function of p for 6 He for the repulsive and attractive neutron- 4 He potentials 
described in the text. The lowest A for the attractive potential originating from zero and corresponding to the Pauli forbidden 
state is removed. 



4 



IV. EFFECTIVE THREE-BODY ADIABATIC POTENTIALS 



e He. The effective radial potential in eq.(|3j) is the crucial quantity which in turn mainly is determined by the A- 
spectrum. The two prescriptions for dealing with the Pauli principle lead to different spectra as seen in fig. for 6 He. 
The purely repulsive neutron-core s-wave potential (solid curves) results in the steep increase of the lowest A-value 
at small p. The increase is quickly interrupted by two (avoided) crossings involving the next two A-values, which in 
contrast are moving steeply down due to the strong p-wave attraction leading to a low lying p-resonance in He (see 
the parameters given in the previous section). They therefore must contain a substantial amount of p-wave and the 
related components must be dominating in the wave function which accordingly also has about 87% p-wave and 13% 
s-wave in the relative neutron-core system. The curves eventually return back recovering the hyperspherical spectrum 
for p = oo. This occurs without (avoided) level crossings. The interaction between crossing levels is vanishing or small 
indicating different symmetries of these levels. The related almost preserved quantum numbers are the neutron-core 
relative angular momenta l x = and 1 corresponding to states originating from zero and from 12. The two levels 
originating from 32 are apparently at small p dominated respectively by s- and p-wave components. At larger p the 
highest A bends over due to p-wave admixture and the related strong attraction. 

The A-spectrum for the other potential with sufficient attraction in the neutron-core s-wave to support one bound 
s-state is also shown as the dashed curves in fig. |l|. The lowest A-value is not shown in the figure. It is mainly s- waves 
in the neutron-core subsystem and decreases at small p and diverges parabolically for p — > oo as clear signals of strong 
attraction and a resulting bound state. All states built on this (not plotted) level are Pauli forbidden and the level is 
therefore excluded in the following calculations. The two next A-values on the other hand remain almost unchanged 
(apart from the avoided crossings) indicating a dominating p-wave content (about 93%) in the relative neutron-core 
subsystem. This behavior is consistent with an almost unchanged energy of the lowest lying p-resonance in the two- 
body subsystem. As the lowest lying levels dominate the wave function the behavior of these effective potentials also 
guarantee an almost identical wavefunction in the two approaches. Also the two highest A-values remain essentially 
unchanged. The third A-value start out as an s-wave as seen from the strong decrease. At larger distances it is then 
replacing the s-wave originating from zero (crossing two times) from the purely repulsive potential. 
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FIG. 2. Angular eigenvalue spectra A n as function of p for n Li for the shallow and deep neutron- 9 Li potentials described in 
the text. The lowest A for the deep potential originating from zero and corresponding to the Pauli forbidden state is removed. 
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ll Li. For n Li the low lying s-state of 10 Li results in a structure where the dominating neutron-core relative con- 
figuration is an s-state. This is rather different from 6 He where the predominant neutron-core relative configurations 
were p-states. The A-spectra for the two potentials corresponding to the different prescriptions are shown in fig. |^. 
The shallow attractive potential without a bound s-state (solid curves) results in a rather slowly changing lowest 
A-value which contains essentially only neutron-core s-waves. This would by far carry the largest probability in the 
resulting wave function. The two highest A-values are also smooth functions of p while the second level steeply 
decreases at small p indicating attraction in the corresponding partial wave. 

The other potential, which has a Pauli forbidden bound neutron-core s-state, has much stronger s-wave attraction 
and a A-spectrum, see dashed curves fig. ^[ with a fast moving level originating from 32. This level responds to the 
large s-wave attraction and crosses quickly the two levels originating from 12 before replacing the lowest level from 
the shallow potential. This level will dominate the configuration of the radial wave function. The lowest A-value 
originating from zero is not plotted and also omitted in the following calculations. It decreases at small p and diverges 
parabolically at large p. The four remaining levels all have similarly behaving counter parts in the other spectrum. 
Again we see the (avoided) crossings indicating almost conserved quantum numbers on the levels. As before they can 
be traced back to the neutron-core relative angular momentum l x — and 1. 
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FIG. 3. The radial wave functions for 6 He (solid curve) and 11 Li (dashed curve) corresponding to the lowest dominating 
A-values for the potentials with Pauli forbidden bound states. 

Radial functions. The radial wave function corresponding to the lowest A is shown in fig. || for the two nuclei. 
This component is responsible for more than 90% of the wave function. Only the result of one of the prescriptions 
is shown for each case, since the curves would be difficult to distinguish especially at larger distances where they are 
very close. This is closely related to the fact that the binding energies and radii are the same in both prescriptions. 
The wave function is spatially more extended for 11 Li than for 6 He, again reflecting the difference in binding energy. 



V. MOMENTUM DISTRIBUTIONS 



The momentum distributions after fragmentation of halo nuclei is a very direct way of gaining detailed information 
about the relative wave function of these nuclei. It was originally expected simply to provide the Fourier transform, 
but more complicated and detailed analyses are needed. Such computations are fortunately available at this moment 
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and we shall exploit this avenue to test our prescriptions. 

The sudden approximation. We consider a process where a high-energy three-body halo projectile instantaneously 
looses one of the particles without disturbing the remaining two. We also assume a light target and we shall therefore 
neglect the Coulomb dissociation process, which then only contributes marginally. 

We work in the center of mass system of the three-body projectile and denote by kj, and k x the total and relative 
momentum of the two remaining particles in the final state. The transition matrix of the reaction in the sudden 
approximation is then given by 

.\/^ , T ;k ; ,.k,.i cx (e <k ^x %ff ^ s ^(k x ,x)|* JM ) , (7) 

where ^ jm is the three-body wave function and w SxC73 . is the final state distorted two-body wave function pd| ] 
corresponding to the two remaining particles with the distance x, total spin and projection equal to s x and o~ x . 
The distance between the center of mass of the two-body system and the removed particle is y and Xs a is the spin 
wave function of the third particle where s y and o~ y are the related total spin and projection. 

The cross section or momentum distribution is now obtained by squaring the transition matrix and subsequently 
averaging over initial states and summing over final states: 

sr^r«E E \Mi M m (^K)\ 2 - (8) 

X " M s x a^a y 

Using the momentum p(= ak^ + bk y ) of one of the particles relative to the center of mass of the projectile as the 

e' 

d 6 a 1 d 6 a 



variable instead of h x , we obtain the relation 



x 



E \Ml M m (K,K)\ 2 , (9) 



dndk. v a 3 dk x dk v 

where a 3 arises from the Jacobi determinant for the transformation. The differential cross section in eq. (§) should be 
integrated over all unobserved variables, i.e. k a and some of the components of p. Note that we have not specified 
any coordinate system, and the axis x, y, and z are therefore completely arbitrary. Thus, in our approximation the 
longitudinal and transverse momentum distributions are identical. 

After neutron removal fragmentation reactions is usual to define the invariant mass E COIC + n as 

£coro+„ = ((^corc + E n f + C 2 ( Pcoro + Pn ) 2 ) 1/2 - (M core + M n ) C 2 , (10) 

where E corre ,m P CO ro m an d AT corCin denote the energy, three-momentum, and rest mass of the core and the neutron, 
respectively. 

Computing the invariant mass in the frame of the two-body system after the fragmentation (p corc + p n = 0) we 
have 

2 2 2 2 _ 1 M COIC M n 2 , _ 2 2 4 _ 2 - . 

Pcoro C — Pn C — ~Tf , ji# K x C ~ ^corc _ JU corc c ~ Al 1Kt n C > k 1L ) 

m M COIC + M D 



and 



dE cmc + dE n = C ° 1C " 7TT^-rTTr k x dk x (12) 
E cole E n m(M coro + M n ) 



The invariant mass spectrum is then defined 

da E COIe E n m(M cole + M n ) 1 da 



where da/dk x is obtained from eq.(pft after integrating over the unobserved quantities. 



(13) 
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FIG. 4. Longitudinal core momentum distributions for a 6 He neutron removal reaction computed for both the purely repulsive 
potential (solid curve) and the attractive potential with one bound s-state (dashed curve) defined in fig. [j]. The broadest 
distributions correspond to the neg lect of final state interaction. The experimental points for fast 6 He (400 MeV/u) colliding 
with a carbon target are from pit], and correspond to the transverse core momentum distribution. The core momentum is 
referred to the center of mass system of the three-body projectile. 
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FIG. 5. Longitudinal neutron momentum distributions for a 6 He neutron removal reaction computed for both the repulsive 
potential (solid curve) and the attractive potential with one bound s-state (dashed curve) defined in the text. The broadest 
distributions correspond to the neglect of final state interaction. The experimental points for fast 6 He (800 MeV/u) colliding 
with a carbon target are from |^|, and correspond to the transverse neutron momentum distribution. The neutron momentum 
is referred to the center of mass system of the three-body projectile. 
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FIG. 6. Radial neutron momentum distributions for a 6 He neutron removal reaction computed for both the repulsive potential 
(solid curve) and the attractive potential with one bound s-state (dashed curve) defined in the text. The broadest distributions 
correspond to the neglect of final state interaction. The neutron momentum is referred to the center of mass system of the 
three-body projectile. 

6 He fragmentation. In fig. ^ the computed longitudinal core momentum distribution after neutron removal in a 
fast nuclear break-up reaction of 6 He on a light target is shown. The results of the two prescriptions are exhibited both 
with and without inclusion of the final state interaction. As expected the effect of the final state interaction is visible 
but not extremely important even with this fairly light core. The results of the two prescriptions are very close when 
the final state interaction is included. When it is neglected the intrinsic differences in the wave functions obtained 
with the two different prescriptions to account for the Pauli principle produce a clear distinction. The difference 
is in the additional node of the forbidden state prescription compared to the reduced attraction prescription. The 
experimental data in fig. ^ are taken from pl[ , and correspond to the transverse core momentum distribution after fast 
6 Hc fragmentation on a carbon target. The present model computations do not distinguish between directions, but 
the results are expected to be more appropriate for the longitudinal directions, which unfortunately are unavailable 
at this moment. However, the agreement with the measured results is still encouraging for several reasons. First, we 
have basically no free parameters. Second, the experimental transverse momentum distributions are expected to be 
broader than the longitudinal ones p2[ . Third, no broadening due to neglected effects and experimental resolution is 
included in the computation. 

In figs. H and ^ we show corresponding computed longitudinal and radial (p 1 - = (p^. + Py) 1 ^ 2 ) neutron momentum 
distributions for a neutron removal 6 He fragmentation reaction. Again the two prescriptions give very similar results 
when the final state interaction is included. In this case we obtain a very large influence of the final state interaction 
which reduce the full width at half maximum by a factor 2-3. Inclusion of the final state interaction is therefore 
necessary to obtain the observed agreement with the measured values. We want here to emphasize that the computed 
results are found in a consistent model where the same two-body potential is responsible for both the three-body 
structure of the initial halo state and the final state interaction after break-up. The experimental data in fig. 
are obtained from p^ j, and correspond to the transverse neutron momentum distribution after 6 He fragmentation 
on a carbon target at 800 MeV/u. Again the computed curve is expected to compare more favourably with the 
experimental longitudinal neutron momentum distribution. However, these data are not presently available. 
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FIG. 7. Longitudinal neutron momentum distributions for a 6 He core break-up reaction computed for both the repulsive 
potential (solid curve) and the attractive potential with one bound s-state (dashed curve) defined in the text. The broadest 
distributions correspond to the neglect of final state interaction. The neutron momentum is referred to the center of mass 
system of the three-body projectile. 



He+Target^n+n+X (X/He+Y) 




p x (MeV/c) 

FIG. 8. Radial neutron momentum distributions for a 6 He core break-up reaction computed for both the repulsive potential 
(solid curve) and the attractive potential with one bound s-state (dashed curve) defined in the text. The broadest distributions 
correspond to the neglect of final state interaction. The neutron momentum is referred to the center of mass system of the 
three-body projectile. 
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Although measured values are not available we show in figs. [7] and |^ the computed longitudinal and radial neutron 
momentum distributions for core break-up reactions of 6 He. In this case the core is violently removed, and the final 
neutron-neutron interaction is relevant. The two prescriptions again lead to similar results. The final state interaction 
can not be neglected due to the low lying virtual s-state in the neutron-neutron system. 

Key quantities for a number of different cases are given in table I, where the 6 He binding energy and radius are 
kept at the measured values. For the case of the repulsive s-wave potential we have first considered two cases (first 
two rows of table I) where both the pi/ 2 and the p 3 / 2 -resonance energies lie within the experimental values. In these 
two cases the computed scattering length for the s-wave is smaller than the experimental value of —3.07. Increasing 
the absolute size of the computed s-wave scattering length (keeping fixed the energy of the lowest 2?3/2-resonance) 
reduces the energy of the pi/2-resonance (case 3 in the table). When the scattering length is equal to the experimental 
value the energy of the second p-resonance is below the experimental energy. Cases 4 to 6 in table I are analogous to 
cases 1 to 3, but with an s-wave repulsive potential giving a clearly worse agreement with the experimental s-wave 
scattering length. However the changes in the full width at half maximum of the momentum distributions, the p-wave 
content and root mean square radius, are very small, reflecting that the properties the s-wave interaction do not play 
an essential role in the structure of 6 He. 

For an attractive s-potential with one bound state (two last rows of table I) it is possible to reproduce both the 
experimental s-wave scattering length and the two p-wave resonance energies. For the two cases 7 and 8 shown in 
table I the variation of the energy of the Vi/i -resonance only produces a change in the 6 He binding energy of about 70 
keV, but in both cases consistent with a computed value of 1.0 MeV. The cases 1 and 7 correspond to the parameters 
given in section 3, and they are the two cases considered for 6 He in figs. |J, ||, [| 0, and ||. 

11 Li fragmentation. The interaction parameters for the schematic n Li computation, where the core spin is assumed 
to be zero, only provide the s- and pi/2-strengths for fine tuning. We already placed the p3/2~state at a very high 
energy and thereby removed it from the active space. In contrast to 6 He the dominating component is now an s-wave 
in the neutron-core subsystem. In table II we show the relevant quantities for three different computations. The 
first row corresponds to a computation with a shallow s-wave potential without bound states while the second and 
third rows correspond to computations with deep s-wave potentials with one bound state. The effective range and 
scattering length are the same in the three cases and the s-wave interaction places the lowest virtual s-state at 200 
keV. Once the s-wave potential has been fixed the energy of the Vi/2 -resonance must then be around 1.7 MeV for 
the shallow potential case in order to fit the 11 Li binding energy. When the deep potential is used we consider two 
situations. First we keep the energy of the p^/ 2 -resonance unchanged. Then the pi/2-energy is reduced to 0.8 MeV 
in order to recover the 11 Li binding energy. In the second case we force the jq^-resonance energy to take the same 
value as for the shallow potential. The correct n Li binding energy is in this case obtained by reducing the energy of 
the p3/2-resonance. The cases 1 and 2 correspond to the parameters given in section 3 and they will be the two cases 
considered for 11 Li in the next figures (more elaborated and detailed computations for n Li fragmentation reactions 
may be found elsewhere, see |IT|). 



case 


E(s 1/2 ) 


a 


E( Pl/2 ) 


r(Pi/ 2 ) 


r c 


r„ 


p-content 


R 




(MeV) 


(fm) 


(MeV) 


(MeV) 


(MeV/c) 


(MeV/c) 


(%) 


(fm) 


1 


5.27 


-2.58 


1.94 


4.02 


101.6 


72.6 


87.7 


2.50 


2 


5.13 


-2.21 


2.55 


7.38 


100.5 


74.8 


86.2 


2.50 


3 


5.09 


-3.07 


1.28 


1.77 


103.0 


68.1 


89.4 


2.50 


4 


10.6 


-2.04 


1.94 


4.02 


98.8 


72.9 


86.7 


2.53 


5 


10.8 


-1.80 


2.55 


7.38 


98.3 


74.5 


85.5 


2.53 


6 


9.3 


-2.51 


1.10 


1.34 


99.2 


67.1 


88.5 


2.53 


7 


5.08 


-3.07 


1.94 


4.02 


107.6 


67.4 


92.8 


2.45 


8 


5.08 


-3.07 


2.55 


7.38 


106.2 


68.1 


92.7 


2.49 



TABLE I. Key quantities for 6 He corresponding to various neutron-core interactions which are chosen as gaussians in each 
partial wave. The neutron-neutron interaction is given in Section 3. For the neutron- 4 He subsystem we give energies of the 
virtual s^-state, s-wave scattering lengths, energies and widths of the p^-resonance. The p 3 / 2 -resonance energy and width 
are in all cases 0.77 MeV and 0.73 MeV, respectively. The full width at half maximum for core and neutron momentum 
distributions is denoted T c and T n . We also give the probability for finding the neutron- 4 He subsystem in a p-wave in the 
three-body wave function of 6 He. The remaining probablity is found in s- waves. The total binding energy of the three-body 
system is in all cases equal to about 1.0 MeV and the root mean square radius is given in the last column. Cases 1-6 and 7-8 
correspond to potentials without bound and with one bound s-state, respectively. 
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Looking at the root mean square radius we observe that it only varies between about 3.30 fm and 3.40 fm indepen- 
dently of the prescription and mainly sensitive to the amount of s-statc in the neutron-core subsystem, see also [[l5) . 
The deep potential gives almost the same radius as the shallow potential. This is in contradiction with the large radii 
(> 3.50 fm) obtained in similar Faddeev calculations, where the lowest s-state is projected out H]. These surprisingly 
large radii are hard to reconcile with the general asymptotic relation between binding energy and radius for weakly 
bound halo nuclei Incidentally a similarly large radius (3.55 ± 0.10 fm) is also obtained in recent Glauber 

analyses of the measured reaction cross section p5|. The reliability of such analyses still remains to be investigated. 




FIG. 9. Longitudinal core momentum distributions for a n Li neutron removal reaction computed for both the shallow 
potential (solid curve) and the deep potential with one bound s-state (dashed curve) defined in the text. The broadest 
distributions correspond to the neglect of final state interaction. The experimental points for fast n Li (468 MeV/u and 648 
MeV/u) colliding with a aluminium target are from [^6|. The core momentum is referred to the center of mass system of the 
three-body projectile. 



case 


E{sy 2 ) 


a 


S(pi /a ) 


r(pi/a) 


r c 


r„ 


p-content 


R 




(MeV) 


(fm) 


(MeV) 


(MeV) 


(MeV/c) 


(MeV/c) 


(%) 


(fm) 


1 


0.200 


8.738 


1.68 


4.16 


57.6 


34.8 


17.3 


3.34 


2 


0.200 


8.738 


0.77 


0.89 


56.1 


38.2 


21.0 


3.35 


3 


0.200 


8.738 


1.68 


4.16 


55.6 


35.6 


15.6 


3.37 



TABLE II. Key quantities for lx Li corresponding to various neutron-core interactions which are chosen as gaussians in 
each partial wave. The neutron-neutron interaction is given in section 3. For the neutron- 9 Li subsystem we give energies 
of the virtual s^-state, s-wave scattering lengths, energies and widths of the p\ / 2 -resonance. The width and energy of the 
p 3 / 2 -resonance is unphysical and not included. The full width at half maximum for core and neutron momentum distributions 
is denoted F c and L n . We also give the probability for finding the neutron- 10 Li subsystem in a p-wave in the three-body wave 
function of 1 Li. The remaining probability is found in s-waves. The total binding energy of the three-body system is in all 
cases equal to about 0.30 MeV and the root mean square radius is given in the last column. Cases 2 and 3 correspond to 
potentials with one bound s-state. 
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FIG. 10. Longitudinal neutron momentum distributions for a n Li neutron removal reaction computed for both the shallow 
potential (solid curve) and the deep potential with one bound s-state (dashed curve) defined in the text. The broadest 
distributions correspond to the neglect of final state interaction. The neutron momentum is referred to the center of mass 
system of the three-body projectile. 
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FIG. 11. Radial neutron momentum distributions for a n Li neutron removal reaction computed for both the shallow potential 
(solid curve) and the deep potential with one bound s-state (dashed curve) defined in the text. The broadest distributions 
correspond to the neglect of final state interaction. The experimental points for fast 11 Li (280 MeV/u) colliding with a carbon 
target are from pSI. The neutron momentum is referred to the center of mass system of the three-body projectile. 
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In fig. ^| we show the computed longitudinal core momentum distributions after neutron removal in a fast nuclear 
break-up reaction of n Li on a light target. We assume that the 9 Li-core has spin zero. The model is then without the 
proper spin couplings and the spin splitting of the neutron-core relative states jig ]. Although simplified the model 
is still not far from being realistic. The two prescriptions again give very similar results and the effect of the final 
state interaction is now a little smaller than for 6 He due to the larger mass of the core nucleus. The experimental 
distributions are essentially reproduced, and correspond to a 11 Li fragmentation reaction on Al at 468 MeV/u and 
648 MeV/u_p|]- 

In figs. |Hpmd we show computed longitudinal and radial (jp 1 - — (p 2 + Py) 1 ^ 2 ) neutron momentum distributions 
for a neutron removal reaction of n Li. The cases shown in the figure correspond to rows 1 and 2 of table II. The 
difference in energy of the p±/2 resonance in the neutron-core subsystem produces the difference in the momentum 
distribution when the final state interaction is included. In fact the momentum distributions obtained in the third 
case of table II are almost indistinguishable from those obtained in the first case (the energy of the P1/2 resonance is 
the same in these two cases) . The final state interaction is then important due to the presence of low lying virtual 
s-states and p-resonances, and very sensitive to the energies of these low lying states. 




FIG. 12. Longitudinal neutron momentum distributions for a 11 Li core break-up reaction computed for both the shallow 
potential (solid curve) and the deep potential with one bound s-state (dashed curve) defined in the text. The broadest 
distributions correspond to the neglect of final state interaction. The neutron momentum is referred to the center of mass 
system of the three-body projectile. 

For completeness we show in figs. |l2]and[l3| the longitudinal and radial neutron momentum distribution from a n Li 
core break-up reaction. We again observe an important effect produced by the inclusion of the final neutron-neutron 
interaction (the neutron-neutron interaction has a low lying virtual s-state), and an almost identical result for both 
the shallow and deep neutron- 9 Li s-wave potential. In figs. |ll| and the experimental data are fairly well reproduced, 
and correspond to a 11 Li fragmentation reaction at 280 MeV/u on a carbon target p9|j27[ | 

Invariant mass spectrum. Computing the invariant mass -E core +n (eq. (|To|) ) in the frame of the two-body system in 
the final state (p corc + p n = 0) we can interpret B corc+n as the kinetic energy of the neutron-core system in the final 
state: 

+ 2/i 2fi 2m V ' 
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where fi and m are the reduced mass of the two-body system and the arbitrary normalization mass, respectively. The 
energy of a resonance is defined as the k\j1m value for which the cross section has a maximum. As a consequence 
the invariant mass spectrum ( |l3| ) will show a peak at the E COTe+n energy equal to a resonance energy. 




FIG. 13. Radial neutron momentum distributions for a 11 Li core break-up reaction computed for both the shallow potential 
(solid curve) and the deep potential with one bound s-state (dashed curve) defined in the text. The broadest distributions 
correspond to the neglect of final state interaction. The experimental points for fast n Li (280 MeV/u) colliding with a carbon 
target are from ^7|. The neutron momentum is referred to the center of mass system of the three-body projectile. 

In fig. [l4|we show the invariant mass spectrum (|l^) from a 6 He neutron removal reaction. The cases with repulsive 
s-potential (solid line) and attractive s-potential with a Pauli forbidden state (short-dashed line) are shown. The 
parameters correspond to cases 1 and 7 in table I, respectively. As seen from the figure, both curves are very similar, 
and show a peak around £/ C ore+n = 0.8 MeV, that corresponds to the energy of the p 3 / 2 resonance. Another p- 
resonance is present at 2 MeV in 5 He. However its width is much larger than before, and therefore its effect on the 
spectrum is much more spreaded out. To illustrate how the invariant mass spectrum is sensitive to the resonance 
energies we also show the case where the neutron- 4 He potential introduced in ref. Q is used (long-dashed line). The 
suggested increase by 3% of the range of the neutron- 4 He potential makes the energy of the p 3 / 2 resonance too small 
(around 0.3 MeV), giving rise to the pronounced peak in the invariant mass spectrum at that energy. The different 
height of the curves in fig. [u] comes from the fact that all the three curves are normalized to 1. 

In fig. ^ we show the same spectrum as in fig. ^ from a neutron removal 11 Li reaction. The cases of the shallow 
and deep potentials are shown (cases 1 and 2 in table II). The difference between these two cases comes from the 
different energy of the lowest p-resonance. 

For the deep potential (dashed line) the pi/ 2 resonance is at 0.77 MeV, creating a clear shoulder in the curve at 
that energy. Note that the contribution to the wave function from the p-wave (21% probability) is given directly by 
the schematic model. It is determined by the energies of the virtual s\/ 2 state and the p\/ 2 resonance together with 
the requirement of fitting the binding energy and radius of n Li. Therefore the shoulder in the distribution is less 
pronounced than in calculations where the p-wave content is chosen to be higher, for instance in psfl , where a rather 
arbitrary 50% contribution for the p-waves is chosen. For the shallow potential (solid curve) the large width of the 
resonance at 1.7 MeV broadens the corresponding peak at that energy and it disappears completely. For both the 
deep and the shallow potential the lowest peak comes from the low lying virtual s-state in 10 Li. However in this case 
the position of the peak is not directly related to the energy of the virtual state. A virtual state produces an increase 
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FIG. 14. The invariant mass spectrum for a 6 He neutron removal reaction computed for both the repulsive potential (solid 
curve) and the attractive potential with one bound s-state (short-dashed curve) defined in the text. The case of the potential 
introduced in ref. (with an increase of 3% of the range of the potential) is also shown (long-dashed line) . 
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FIG. 15. The invariant mass spectrum for a 1 Li neutron removal reaction computed for both the shallow potential (solid 
curve) and the deep potential with one bound s-state (short-dashed curve) defined in the text. In the inset the case of the deep 
potential (solid curves) and the contributions from the s-waves (dot-dashed line) and the p-waves (dashed line) are shown. The 
experimental points for fast n Li (280 Mev/u) colliding with a carbon target are from p9| . 
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of the momentum distributions at zero momentum. In fact, the value of the invariant mass spectrum divided by 
\J -Ecore+n is not zero at zero energy, and the lower the virtual s-state the larger the value at the origin. In the inset 
of fig. [l5| we show the deep potential case where we have separated the contributions from the s and the p-waves. It 
is then clear how the s-wave is responsible for the first peak, while the shoulder is produced by the pi / 2 resonance at 
0.77 MeV. The experimental points |2j| correspond to a n Li fragmentation process on carbon. A better agreement 
between the computed curves and the experimental data is possible when the spin of the core is taken into account 
in the description of 11 Li and 10 Li. 

VI. SUMMARY AND CONCLUSIONS 

We study two prescriptions to take the Pauli principle into account in the three-body cluster model where more than 
one cluster contain nucleons. One is to exclude the Pauli forbidden two-body state from the active space available for 
the three-body system, the other is to construct a two-body potential without the Pauli forbidden state but with the 
same low energy properties. 

We solve the Faddeev equations in coordinate space by means of the adiabatic hyperspherical expansion. We 
calculate the angular eigenvalues of the Faddeev equations which are closely related to the effective radial potentials. 
This spectrum is computed for the two different potentials with identical low energy properties and either without or 
with one Pauli forbidden bound s-state. Each bound state gives rise to a specifically behaving angular eigenvalue. The 
technically inexpensive prescription to account for the Pauli principle is then simply to omit the angular eigenvalue 
which corresponds to the Pauli forbidden bound state in the calculations of the three-body radial wave functions. 

We apply these prescriptions to analyses of the high energy fragmentation reactions of two halo nuclei 6 He and 11 Li 
within a three-body neutron+neutron+core model. The potentials of the two prescriptions are adjusted to reproduce 
the same s- and p-wave low energy data. 

For 6 He these data are known experimentally and the interactions are therefore fixed. The potential with forbidden 
state then gives correct binding energy and root mean sqare radius for 6 He. The equivalent potential without a bound 
state slightly underbinds the system. A small reduction of the repulsive core is needed to obtain the correct binding. 
This automatically provides almost the same reasonable root mean square radius. After this small adjustment the 
wave functions for the two nuclei for the two types of potential are compared and found remarkably similar. 

For 11 Li we use the currently accepted experimentla data for the lowest s- and p- levels in the neutron core subsystem 
to adjust the neutron core potentials. With these levels being close to experimental data the deep potentials provide 
reasonable binding energy and root mean square radius. The equivalent potential without a bound state slightly 
over binds the system. After a small attenuation of the position of the neutron core P3/2 level this prescription also 
provides the correct binding and size of the system. 

In general the two prescriptions without any fine tuning provide very close but still distinguishable ground state 
properties of the three-body system. After small adjustments the properties become remarkably similar. 

We then compute neutron and core momentum distributions in nuclear break-up reactions of these halo nuclei. 
We use the sudden approximation and include final state interactions, which are crucial for the neutron distributions 
as the full width at half maximum for the neutrons from 6 He is reduced by a factor of more than 2 from the value 
corresponding to omitting the final state interaction. 

The two prescriptions give almost identical results and the measured momentum distributions are rather accurately 
reproduced. These computations are carried out in a consistent model where the same two-body potential is used 
both for the initial three-body halo structure and for the final state interaction after the break-up process. 

The final state invariant mass spectra for 5 He and 10 Li are computed and their features discussed. A virtual s-state 
shows up as a pronounced peak close to zero energy. This peak arises from the phase espace factor which vanishes at 
zero energy, and its position is unrelated with the energy of the virtual state. Higher angular momenta show up as 
distinct peaks at the resonance energy. These spectra are then sensitive to the continuum structure of the two-body 
system. 

In conclusion the two prescriptions to account for the Pauli principle work remarkably well for the two test examples 
6 He and 11 Li. In both cases the lowest Pauli forbidden s-state is occupied by nucleons in the core nucleus. However, 
the valence neutrons occupy predominantly neutron-core relative p-states and s-states, respectively for 6 He and n Li. 
The wave functions are further tested in connection with our fragmentation model and the measured momentum 
distributions are nicely reproduced. For 6 He these quantities are computed for the first time whereas n Li was 
investigated earlier in a more sophisticated model. Although fairly realistic, we used here only a schematic model, 
because the main purpose was to demonstrate the reliability to account for the Pauli pinciple. In any case the results 
provide additional support both for the prescriptions and for the fragmentation model. 
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